Optimizing local protocols implementing nonlocal quantum gates 
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We present a method of optimizing recently designed protocols for implementing an arbitrary 
nonlocal unitary gate acting on a bipartite system. These protocols use only local operations and 
classical communication with the assistance of entanglement, and are deterministic while also being 
"one-shot" , in that they use only one copy of an entangled resource state. The optimization is in the 
sense of minimizing the amount of entanglement used, and it is often the case that less entanglement 
is needed than with an alternative protocol using two-way teleportation. 

PACS numbers: 03.67.Ac 



I. INTRODUCTION 



Recent years have seen enormous advances in the study of quantum information, including in specific areas 
such as quantum computing [U-IH. O ne particular proposal for a quantum computer attempts to overcome 
the difficulty of protecting a large quantum processor from decoherence by using a multi-processor design 
so that the size of each processor may be kept small. Obviously, such a "distributed" quantum computer 
Q will require the implementation of tasks that involve multiple processors, including information transfer, 
measurements, and other nonlocal operations utilizing unitary gates. One possibility for performing such 
operations is to physically move component quantum systems around, but it may well be advantageous to 
perform these operations without bringing parts of different processors together, but rather by utilizing the 
power of entanglement. In such a situation, one is considering spatially separated parties restricted to local 
■ quantum operations on their individual subsystems with communication of classical information amongst 

<3\ \ the parties (LOCC). Although numerous tasks exist that remain possible under the LOCC paradigm when 

entanglement is not available — an important example involves distinguishing between two orthogonal 
quantum states Q — there are also many important cases where global operations are required. However, 
even the latter tasks become possible with LOCC when there is enough entanglement available to the parties. 
This is one of the main reasons why entanglement is now recognized to be a valuable resource, especially 
when considering that it is difficult to create and maintain. An important question, then, is this: For a given 
nonlocal task, how much entanglement is necessary? 

Due to the importance of such questions, much attention has been focused on the design and optimization 
of protocols utilizing entanglement to perform nonlocal tasks Q-QjJ • Such protocols fall into a variety of 
classes, including those that utilize an asymptotically large number of inputs in order to obtain a fraction 
of outputs having vanishingly small error [71 Il5|. In contrast, deterministic protocols require perfect success 
with each individual input. In this case, one may consider that the given task will be repeated many times 
with the goal to optimize entanglement resources in terms of the average amount used over all repetitions 
[1, [OJ. On the other hand, one may also consider one-shot protocols [lQ, El, EH El|, m which the task is 
only performed once. Then, one wishes to know the smallest amount of entanglement with which the task 
may be accomplished in a single implementation. While these latter protocols will generally require the 
most entanglement, they do offer the advantage of certainty (which may at times be needed) and as such a 
savings in other resources, including the amount of time needed for the task. 

Here, I consider the one-shot deterministic implementation of nonlocal unitaries acting on bipartite system 
Ha <8> W-B using only LOCC with the assistance of entanglement. We have recently succeeded in creating 
efficient protocols, with the most general of these having the ability to implement any bipartite unitary [161 ]. 
The crucial remaining question is to determine the level of efficiency, in terms of the amount of entanglement 
that is required, that can be attained with these protocols. One way to measure efficiency is by comparison 
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to a "standard" protocol, in which the state of one subsystem, say Ha of dimension Da (let us assume 
Da < Db), is teleported to the laboratory of the other party, who performs the desired nonlocal unitary 
on the combined systems before teleporting the final state of the first party's subsystem back to where it 
belongs. This standard protocol requires two maximally entangled states (MES), each of rank Da, so uses 
an amount of entanglement equal to 2log Da- We thus wish to know under what circumstances, and in 
particular for which nonlocal unitaries. our protocol succeeds while using less entanglement than this. This 
question was essentially answered in [lfjj, where a method was provided that parameterizes all bipartite 
unitaries that can be implemented by our protocol using a given amount of entanglement. However, it was 
not known how to determine the minimum amount needed for a particular, given unitary, a question of 
crucial importance since circumstances will often dictate which specific unitary is needed. The purpose of 
the present paper is to provide an answer to this question. 

The protocol of [la ] for implementing the nonlocal unitary 11 is based on an expansion of the form 

u = J2u(f)(sw(f), (1) 

fee 

where matrices U(f), which act on Ha, make up a projective unitary representation of a finite group G. 
The W(f) are arbitrary matrices, apart from the requirement that U be unitary, and act on Hb- To 
deterministically implement this unitary, our protocol uses — and in fact when G is the smallest group 
allowing an expansion of U in the form ([T]), requires [TtJ — an amount of entanglement equal to log 
where \G\ is the order of the group G. Given 14, it is always possible to arrive at such an expansion by 
taking the U(f) to be the generalized Pauli operators of dimension Da (for a definition, see Eq. (11) of 
[ltl). However, since the order of this group is D A , this particular expansion uses the same amount of 
entanglement as does the standard protocol, so it will not generally be optimally efficient. The reason this 
group may always be used is that the given representation in terms of generalized Pauli matrices forms a 
basis of the space of operators acting on Ha- This will obviously not be true if one chooses a smaller group, 
and for this reason an arbitrary group will not generally allow an expansion of a particular, chosen U. So 
given an arbitrary 14, we seek an "optimal" expansion; that is, an expansion in terms of the smallest group 
possible. In the remainder of this paper, we will see how such an expansion may be found. 

We begin in the following section with a description of the protocol of [16|. Then, in Section Hill we show 
how to find the optimal expansion for an arbitrary 14, which involves choosing the smallest group G. We 
conclude with a discussion of our results in Section IIVI 



II. THE PROTOCOL 



The protocol we consider in this paper has been extensively discussed in [16|. Here, we give only a brief 
description. It is based on the use of a finite group G, which may or may not be Abelian, having elements 
/, g, etc., with e the identity. Group multiplication is indicated by fg. The unitary operators U(f) of ([1]) 
form a projective representation of G in the sense that 

U(f)U(g)=fi(f,g)U(fg) (2) 

for all / and g in G. Here the ^{f,g) are nonzero complex numbers constituting a factor system; in our case 
they have unit magnitude because the U(f) are unitary. Without loss of generality, one can always take 
A t (/j<?) = 1 whenever / or g is equal to e or when g = / , and we will follow this convention throughout. 
In addition, associativity of group multiplication requires that 

Kf>9)v(f9,h) = (J,{f,gh)n(g,h). (3) 

A circuit diagram for the protocol is given in Figure [T] Alice and Bob share an MES, 

l$) = T^E |/>«|/>6, (4) 

1 1 fee 

of Schmidt rank \G\ on systems a, b described by Hilbert spaces H a ,Hb- Alice begins the protocol by 
implementing a controlled set of unitaries, where system a in standard basis state \ f) a controls the unitary 



3 



U(f), which operates on Ha- She then measures a in a basis unbiased to the standard basis and defined by 
unitary operator F, where by unbiased we mean that when represented in the standard basis, F must have 
all entries of magnitude | C | 1 / 2 . Next, Alice tells Bob the outcome h of that measurement, and he follows 
by performing diagonal unitary Z(h) on b, which has the effect of removing phase factors introduced by 
Alice's measurement. He then performs unitary M on bB, which introduces the operators W(f) of (J| in a 
way that is correlated with the Z7(/)'s of Alice's earlier operation, where these correlations are made possible 
by the initial entanglement between a and b. Bob then measures b in the standard basis and tells Alice his 
outcome, g. Alice completes the protocol with the "correction" U(gy on A, which adjusts the correlation 
between W's and U's so that the result is always equal to U (that is, so that W(f) is always tensored with 
[/(/), rather than with U(f') for some /' ^ /). 
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FIG. 1: Circuit diagram illustrating local implementation of nonlocal unitary U = ~}2f eG U(f) <g> W(f), where the 
set of unitaries {U(f)} forms a projective representation of a group. 

The matrix M is given by 

M = R(f) ® W(f) (5) 
fee 

with 

R{f) = J2K9J)\9)(9f\, (6) 

and these R(f) form a projective regular representation of G with factor system fi(f,g). It can be shown 
fioi ) for the circuit of Figure [1] that the unitarity of M implies that of U. By using the basis {\g)} of T-Lb 
one can view M as a matrix with blocks, where the (g\M\f) block is equal to W{g _1 f) multiplied by a 
phase [i(g, g -1 /). This structure is responsible for the effectiveness of the correction U (gy in completing the 
protocol so that it successfully implements U for every measurement outcome g by Bob (and h by Alice). 
The reader is encouraged to consult [l6| for more details and further discussion. 



III. FINDING THE OPTIMAL EXPANSION 



In this section, we show how to find an optimal expansion of the form ([T]) for an arbitrary unitary U; 
that is, an expansion in terms of a unitary representation U (/) of the smallest possible finite group G. 
Recalling that the point of finding the smallest group is to minimize the amount of entanglement required 
to implement U using the protocol described in the previous section, we will consider local unitaries to be 
"free" . That is, instead of {[T|), we wish to consider the more general case of an expansion in the form, 

U={V' A ® V' B ) r £ U(f) ® W(f)J (V A 9 V B ), (7) 
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with unitary Va{b) an d Va(b)- ^ or a gi ven group G, there will be unitaries that can be expanded in this way 
but not in the form (JXJ) . The local unitaries acting on Hb simply lead to a replacement, W(f) — > V b W(J)Vb : 
so these may be absorbed into the definition of W(f) and effectively ignored. In addition, 

v A u(f)v A = v A v A (ylu{f)v A ) , (8) 

and since V A U(f)V A is a representation of G equivalent to the representation by £/(/), we may replace 
^(/) Thus, defining y = V A V A , we will include all possibilities by looking for an expansion 

of U in the form, 

u = Y,[yu{f)\®w{f). (9) 

fee 

[Note that, by the above arguments, it would be equally valid to place the V to the right of the U (/) in this 
expression.] 

Therefore, given an arbitrary bipartite unitary 1A, we now see that our task of finding an optimal expansion 
requires the identification of: (1) the smallest group, G; (2) the projective representation, U(f); (3) the local 
unitary, V; and finally, (4) the set of operators W(f) acting on Hb- It should be clear that if G is unknown, 
there is little hope of identifying the other items in this list. By first developing a method of finding G, we 
will see that the others follow without too much additional difficulty. 

We note that one will also want to look at expanding in terms of a group on the % b side, as this may lead 
to a smaller group than the best expansion on the T-i A side even when D A < Db- Of course, if one finds an 
expansion with a group of order equal to the operator Schmidt rank of U, which is the number of terms in 
its Schmidt expansion (see Eq. (|10p . below), then one is done since (by definition) no expansion can have a 
smaller number of terms than this. 

Ill A. The optimal expansion of U 

We begin with the observation that every group is characterized by the dimensions of its irreducible 
representations (henceforth referred to as irreps). Our approach will then be to identify the set of irrep 
dimensions the group must contain, and then to find the smallest group G that contains those dimensions. We 
know from group representation theory that the set {U(f)} in (0) can be reduced — that is, simultaneously 
brought to block-diagonal form — by a unitary similarity transformation, and that the resulting blocks will 
each be associated with one of the irreps of G. Let us now see how we may identify the necessary irreps. 

In order to make progress, we first need a way of expressing U. Let us use its operator Schmidt decompo- 
sition, 

U = Y,A j ®B j , (10) 

3 

where Tr(SlSj) = Sjk, and Tr^AiAj) = Vj : ^ k. Given U, we can always find the Aj and Bj, so we may 
assume these matrices are known. Equating the expansions in (|9|) and (|10p . multiplying by B' k , and tracing 
over Hb, we find that 

v^A k = m f u(f), (ii) 

feG 

with W kf = Tr(4w(/)). 

If we express (ITTt in the appropriate basis the U (/) will all be block diagonal, and then so must be the 
product, V^Ak, for every k. If V was known, then we could block-diagonalize the set of (known) matrices, 
{yUj, and by relating these blocks to those of the U(f) in (fTl]) . we could identify the size of the irreps 
that must be present in any group that the U (/) could possibly represent. However, it is the Ak that are 
known to us, V and U (/) being what we are trying to find. This means that (|11[) is not directly useful for 
identifying the irreps G must contain. What we seek is a known set of operators that can themselves be 
block-diagonalized, and at the same time expressed as a linear combination of the £7(/)'s. The problem in 
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the previous equation resides in the presence of V, which is unknown, so we will use a little trick to get rid 
of this operator. 

Consider, then, multiplying (jlll) on the left by its Hermitian conjugate for index j, yielding 

A)A k = £ /i(g-\/)W*,}V fe/ C/Qr7), (12) 

where use has been made of ©. The objects on the left-hand side of this equation are known, and if the 
U (/) are block-diagonal, then the A^A^ will be so, as well. In fact, the above equation tells us that the finest 
block- diagonal form of the set {U(f)} cannot be finer than that of the set {Aj^4/c}- If we can find the finest 
block-diagonal form of {Aj^lfc} and then choose the smallest G that contains irreps of the same size as those 
resulting blocks, this will be the smallest G that allows an expansion of the form (|12l) . However, there is one 
potential problem with this approach, indicated by (jTTJ) . This equation tells us that the finest block-diagonal 
form of {U(f)} also cannot be finer than that of the set {W^4fc}, and it is not immediately obvious that the 
just-described method of choosing G will necessarily satisfy this additional constraint. That is to say, it may 
be that there is no V such that the set {V^^,} will have as fine a block-diagonal form as the set {A^Afc}. 
However, the following theorem tells us that for sets of operators of the type we are considering, this is not 
a problem. 

Theorem 1 Given a set of operators {A k } taken from an expansion of a non-singular operator in the form 
Ak <E> Bk, there always exists a unitary V such that the set {V^Ak} has the same finest block diagonal 
form as the set {AjAk}. 

The proof of this theorem, which also provides a method of constructing V, is given in Appendix [A"l 

We therefore need to find the finest block-diagonal form of the AjA^. It turns out that a very efficient 
method has recently been devised that accomplishes this task , providing the sizes of the necessary 

irreps that G must contain, as well as the unitary S that block-diagonalizes these operators. Thus, given U, 
we can find the size of the irreps that G must have in order that an expansion of the form ^ is possible. 
Note that the number of different irreps needed is equal to the number of inequivalent blocks in the AjAk- 
For example, suppose that for every j, fc, the finest block-diagonal form is 

4*-(f <)• (13) 

where the two blocks are of equal size d and R^jJ = TR^T^ for some fixed unitary T. Then these two 
blocks are "equivalent" to each other, and there is only a single irrep needed in G, of size equal to d and 
repeated twice in each U(f) [201 ]. 

An important implication of the proof of theorem [T] is that the unitary S which block-diagonalizes the 
AjAk, and which is given by the method of [HI, will also block-diagonalize the V^Ak- Then, once the 
group G is chosen (see below), we will define the representation U(f) to be block diagonal in the same basis 
as are the V^A^ with the appropriate irreps in the appropriate positions. Then, since the proof of theorem 
Q] provides a construction of V, we will have identified V and the representation matrices U{f). What is left 
is to find the operators W(f), and then we will return to the question of actually choosing G. 

To find W(f), we note that 

U = 5>J® B J=£(£ WjfVUtf) ] ® Bj 
j j V/eG 



= £[W(/)]®Wm (14) 
fee 



G 



where in the first line we used (fTTj) and the last line is just ©, allowing us to identify the choice, 

W(f) = J2^fB,- (15) 

3 

Matrix W may be found as follows. Label the inequivalent blocks of V'Ak by the irrep each is associated 
with, so that by (TlT]) . 



(yU fe ) (A) = ^w fe/ c/W(/), 

fee 



(16) 



with C/< A ) (/) the A*' 1 irrep of G. Define 



,(A) 



(17) 



of irreps [21 



where the sum is over all the k inequivalent irreps belonging to G. By the well-known orthogonality condition 

(18) 



Ct^X ' _ V / n'm' \ J mn 



feG 

we see that this definition of W automatically satisfies (fT6|) and if every irrep of G appears in the represen- 
tation by the U(f), we are finished apart from the actual choice of G. If, on the other hand, any irreps A 
of G are missing in the U{f), V^Ak will have no block (V^Ak)^ associated with these irreps, and for the 
purpose of calculating W from (|17p . one may define the latter quantities in any way one wishes, including 
setting them to zero. This is consistent with the fact that missing irreps correspond to linearly dependent 
representations (l6j . and therefore a certain amount of freedom in the expansion coefficients W(f) (compare 
(fT5f ). Note that given any set of W(/)'s, then by theorem 5 of [l6[ there is always at least one choice of 
coefficients for which the protocol described there will work. Thus, given any expansion in terms of the 
U(f), such as one with the W(f) obtained by setting (V^Ak)^' to zero in (fT7|) for all the missing irreps, a 
method of choosing a workable set (that is, one for which our protocol will work) from the original set of 
W(f)'s in the linearly dependent case is described at the end of Appendix D in 16]. Hence, we have shown 
how to obtain an expansion of any IA in terms of a group, an expansion that is optimal for our protocol once 
the smallest group G is chosen. Let us now turn to the task of finding the smallest group. 



Ill B. Choosing the smallest group 



We are left with the task of identifying the smallest group G that contains a set of irreps of given dimensions, 
{d\}, where A labels the required inequivalent irreps (for convenience, order these as d\ < e?2 < ■ • ■)• Here, I 
will describe how this can be done if confining oneself to ordinary representations; the (difficult, but probably 
often necessary) task of including projective representations will be discussed in Appendix [Bj Now, since the 
order of a group, \G\, is equal to the sum of squares of the dimensions of its inequivalent irreps, it must be 
that \G\ > No = d\ , where here represents a sum over only those inequivalent irreps we have found 
are needed in G. Note also that d\ = X^( e )> that is> the irrep dimensions are equal to the character of the 
identity element in that irrep. So, one approach would be to search through character tables of groups of 
order No, looking at characters of the identity element in the various irreps for these groups. If there exists 
such a group with characters of the identity element matching the desired set {d\}, then we are done and 
that group can be used for an optimal expansion of U. Otherwise, we need to look at groups of higher order 
N > Nq, but recalling that the irrep dimensions are divisors of \G\ (a fact that holds even for projective 
irreps |22[), then given {^a}, this provides a strong constraint on possible values of \G\. So considering this 
constraint we can search through groups of higher and higher order until wc find one containing all the 
required irreps. However, if before finding such a group we reach a value of N — (di + e^) 2 + X)a>2 ^A' 
then one must also look for groups that have irrep dimensions, di + d$, ■ ■ -, along with looking for those 
with di, d^ ■ ■ ■■ The reason is that the two blocks of dimensions di and d2 may be considered as a single 
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block of dimension d± + d®, with this larger block associated with a single irrep. This procedure may be 
continued to higher order groups, combining blocks as one proceeds (of size d\ + d$, for example, and so 
on) until one finds a suitable candidate. However, if one reaches N = D 2 A before finding a suitable group, 
one may conclude that it is not possible using the protocol of Figure [T] to implement U more efficiently than 
by teleporting both directions, unless one can find a suitable projective representation of a smaller group, 
as described in Appendix [Bj (Of course, one will also want to check to see if a smaller group can be found 
when the roles of Alice and Bob are reversed, the group representation matrices U (/) then appearing on the 
Hb side instead of Ha-) 

[Another possible approach is the following. For each group satisfying the constraints mentioned in the 
previous paragraph, starting once again with the smallest, construct the regular representation in its standard 
form (permutation matrices with each row and column having only a single non-zero entry equal to one). 
This can be done in a straightforward way from the multiplication table of the group. Then the method of 
[IH [l9| may be used to decompose this set of matrices into irreps, in the same way as was done above for 
the A^-Ak- Note that it is only necessary to construct the regular representation for a set of generators of 
the group and then to decompose these, since if and only if the generators are in their finest block diagonal 
form, then the whole group will be so as well. This is an efficient method of finding not only the dimensions 
but also the irreps themselves, which will eventually be needed so that we can construct the U(f). However, 
if a more efficient method (perhaps, for example, the one described in the previous paragraph) is able to 
determine the dimensions of the irreps of a group, this would be preferable, as we only need to know the 
irrep matrices of the group we ultimately choose to use. The irrep matrices for those groups that are rejected 
in the process are of no interest to us.] 

These searches may be computationally demanding for large dimension, especially when looking for pro- 
jective representations as described in Appendix iBl but will certainly be tractable if Da is not too large. 
We make no claim that the above methods (or that described in Appendix [B] for projective representations) 
are optimally efficient, and we leave the question of optimizing this search as an open problem. 



Ill C. Starting from expansions other than the Schmidt expansion 

There is one final issue that must be addressed. The arguments of the previous section are based on a very 
particular starting point, the Schmidt expansion of U. We must consider whether using a different starting 
point could lead to an expansion in terms of a smaller group. The answer is no, as the following discussion 
makes clear. 

Suppose we start from a different expansion, 

U = Y,A m ®B m . (19) 

m 

Equating this to (1101) . the Schmidt expansion of hi, multiplying by and tracing over T-Lb, we find 

Ak = ^2PkmA m , (20) 

TO 

with (3 km = Tr(S|B m ), and 

A]A k = J2PkmP* n AtA m . (21) 
mn 

This tells us that the finest block-diagonal form of the set {A^A™} cannot be finer than the finest block- 
diagonal form of the set {A^A^}. We may therefore conclude that no other starting expansion can lead to 
a finer block-diagonal form than when starting from the Schmidt expansion of hi, and therefore the group 
obtained by the arguments of the previous section is in fact the smallest group possible. Indeed, a moment's 
thought leads to the conclusion that any expansion of the form (|19[) with linearly independent operators, 
{Bra}, will lead to the same finest block diagonal form as does the Schmidt expansion, and hence also to 
the smallest group (this could be useful, as it may be easier to obtain some other expansion, such as in 
terms of generalized Pauli operators, than to obtain the Schmidt expansion). The reason is that when the 
B m are linearly independent, one may also obtain A m as a linear combination of the A k , so that the finest 
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block-diagonal form of the set {Atylj.} also cannot be finer than the finest block-diagonal form of the set 
{AjjAm}. The conclusion follows immediately, and it is therefore acceptable to use any starting expansion 
one chooses, as long as the {B m } are a linearly independent set. 

IV. CONCLUSIONS 

We have presented a detailed method for finding an optimal expansion of a nonlocal unitary allowing 
IA to be implemented with the smallest amount of entanglement possible when using the protocol of [16]. 
This method utilizes an expansion of U in terms of a projective representation of a finite group G, where the 
order (size) of G determines in a precise way the amount of entanglement used in the protocol. Therefore, 
finding the smallest group is crucial, a task that is likely to be computationally intensive in many cases. 
We have presented possible ways to find this smallest group but leave as an open problem the question of 
optimizing this search. 

The protocol considered in this paper is just one that has been presented in |l6j, and it is possible to 
identify the others as special cases of this one. If, for example, the U(f) all commute with each other, then 
Li is a controlled-unitary (l6j . which can be written in the form 

W = 5^Pj<8V,-, (22) 

3 

with Vj an arbitrary set of unitaries on T-Lb and the Pj are projectors on Ha satisfying P-j = I a- 
This expansion as a controlled-unitary may be found directly from the expansion in terms of the U(f) by 
simply diagonalizing the latter and identifying the Pj as projectors onto the subspaces corresponding to the 
individual one-dimensional irreps. Note that some of these irreps may be repeated in the U(f) representation, 
and this will lead to higher-rank projectors Pj [l6[ • 

The protocol of [16| for controlled-unitaries may have advantages over the protocol discussed here, as the 
local gates used by Alice and Bob may be easier to implement in practice (the unitary M of Figure [T] is 
replaced by a simpler set of gates). This may also be true of the protocol for what has been called a "double 
unitary" representation in [161 ], with an expansion of the form 

U = <f)U{f) ® V(f), (23) 
fee 

where {U(f)} and {V(f)} are each a projective representation of the same group G and c(/) are complex 
coefficients. If the set of operators W(f), found from (fTS)) . are proportional to unitaries that represent G, 
then we have this type of double unitary representation. However, if the U(f) are linearly dependent, then 
even if the W(f) are not proportional to unitaries, there may be another set that can be used that are. 
Finding the easiest way to determine if this is the case is left for future study. 
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Appendix A: Proof of theorem [T1 

Here I prove that if the are the expansion coefficients on one side for a non-singular bipartite operator, 
A4 = J2j Aj ® Bj , then there always exists a single, fixed unitary V such that the finest block diagonal form 

of the V^Ak is the same as that of the A^A^- The argument provides a direct method of constructing V. 

Given that A^Ak are block diagonal for every j, k (assume, without loss of generality, we work in the 
appropriate basis to bring them to this form), with blocks of size rii, ri2, etc., in that order. Then this means 
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the following: The first n\ columns of Ak are orthogonal to columns n% + l,ni + 2, • • • , Da of every Aj, 
including Ak; columns ni + 1, • ■ • , ni + ri2 of Ak are orthogonal to columns n± + ri2 + 1, • ■ • , -Da of every 
A,- (and also to the first n\ columns, of course); etc. So we can partition all the columns of all the Ak into 
subsets of vectors, where any vector in one set is orthogonal to every vector in any other subset. 

As M. is non-singular, the collection of all these vectors spans Ha (Alice's Hilbert space). If this were 
not the case, then there would exist a state \a) S H a with the property that for any \b) £ Hb, the ket 
|£) = | a) <g> 1 6) lies in the nullspace of M\ contradicting the assumed non-singularity of M. Furthermore, 
we have the stronger statement, 

Lemma 2 For any subset of the integers 1 to Da, call this subset S having cardinality n, collect all the 
vectors forming columns mi,m,2, • • • , m n , with rrik 6 S, from all the Aj taken from an expansion of a non- 
singular (bipartite ) matrix M. — ^2j Aj ® Bj . Then, this collection of vectors spans a space of dimension at 
least n. 

Proof of lemma: Defining states |1), ■ • • , \D a) as the standard basis of Ha (the subsystem upon which the 
Aj act), we have for example, that Aj maps state \m\) to the vector defined by column mi of Aj. If the 
statement of the lemma were false, than the collection of vectors {\mk)\l)} (where k = 1, • • • , n; / = 1, • • • , Db] 
and \l) are the standard basis of Hb), which together span a subspace of dimension nDs, would be mapped 
by M. into a subspace of dimension strictly less than nDs- This is a contradiction, since a non-singular 
matrix maps subspaces into subspaces of equal dimension, proving the lemma. H 
Due to this lemma, we see that the partitioning mentioned earlier is such that every subset of vectors spans 
a subspace of dimension equal to the number of columns it is drawn from (n a for the a th subset). No subset 
can span a larger space, since this would mean, considering all the subsets together, having more than Da 
orthogonal vectors on Ha- Choose a basis for each of these subspaces, call these bases B\ for the subset of 
the first n\ columns, B 2 for that of the second ri2 columns, etc. Then we may conclude that every Ak maps 
{|1), • • • , |ni)} into Bi, {\ni + 1), • • • , \ni + n 2 )} into B 2 , and so on. Define V to be a unitary that mirrors 
these mappings. Then, undoes the action of each of the Ak on these subsets of basis states, mapping Bi 
onto {|1), • • • , \ni}}, B 2 onto {|m + 1), • • • , \n x + n 2 )}, etc., and V^A k maps {|1), • • • , |m)} to {|1), • • • , K)}, 
{|ni + 1), • • • , I Hi + ^2)} to {|m + 1), • • • , \ni + n 2 )}, and so on. This means that V^Ak is block diagonal in 
the same form as all the and we are done, including having constructed the unitary V (as the above 

argument indicates, the choice of V is generally not unique). Note also that these arguments imply these 
two sets of operators, {W^4 fe j. anc i {AjAfc}, take on their respective finest block diagonal forms in one and 
the same basis. 

Appendix B: Searching for projective representations 

We wish to extend our search for the smallest group to include projective representations of a group G 
(23l | . If for each admissible N , no group of this order exists that has the desired ordinary irreps, we then 
need to determine if any such group exists that instead has a set of projective irreps matching the required 
dimensions. Note the additional requirement that all irreps must share the same factor system, since if the 
representation by U(f) is decomposed into a direct sum of irrep blocks, then @ will be satisfied if and only 
if it is satisfied block- by-block with a fixed fi(f, g). This point indicates that if d\ = 1 for any A, then it must 
be that the representation is an ordinary one, since one-dimensional projective irreps are always equivalent 
to ordinary irreps. 

Below, we will describe one possible method of searching for a group having the desired irrep dimensions, 
which will utilize the following lemma. 

Lemma 3 A projective irrep of dimension d exists for a group G if and only if an ordinary irrep of that 
dimension exists for a central extension L of G = L/K by some cyclic group K of order r, where r divides 
\G\. Furthermore, all projective irreps of G that share the same factor system {(i(f,g)} may be derived from 
ordinary irreps of a single such extension, L. 

Proof: The "if part follows directly from the well-known relationship between central extensions and pro- 
jective representations 23]. To prove the converse, suppose a projective irrep T of dimension d and factor 
system fx exists. Then there also exists a projective irrep T of this dimension having a normalized factor 
system ll [23]. That is, for all / and g, 



T(f)T(g) = Ll (f,g)T(fg), 



(Bl) 
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where by normalized we mean that fi(f, g) = uJr U " g) with uj r = e 2 ^ r for some r dividing \G\ and n(f,g) 
is an integer. Consider an extension L of G with elements (m, g) 6 L, where g € G, m is an integer lying 
between and r — 1, and group multiplication is defined as 

(l,f)(m,g) = (l + m + n(f,g),fg), (B2) 

with addition mod r. Given our definition of n(f,g) in terms of the normalized factor system /i(/, <?), it is 
not difficult to show that this multiplication satisfies the properties of a group. 
Next, construct from the irrep T(g) the function 

V{m,g)=^T(g), (B3) 

which satisfies 

V{1, f)V(m,g) = u; l r +m T(f)T(g) = J 7 +m+n ^T(fg) = V(l + m + n(f, g), fg), (B4) 

so is clearly an ordinary irrep of L having dimension d, and it is readily seen that L is a central extension 
of G by the cyclic group of order r. Furthermore, we see that g) defines n(f, g) which in turn, along 
with the given group multiplication defined for G, defines the group multiplication for L. Therefore, these 
arguments show that every projective irrep r of G that shares this factor system /i gives rise to an ordinary 
irrep T> of this single extension L, completing the proof. ■ 

One way to determine all the projective irreps of a group G is to find the Schur cover or representation 
group of G [23j . but this is generally a difficult task and in any case, provides much more information than 
is needed here. Let us now describe another possibility, which may be more efficient. We seek a group G 
of order N with projective irreps of specific dimensions {d\}, all sharing a single factor system. By the 
previous lemma, we know that each such irrep exists iff an ordinary irrep of that dimension exists for a 
central extension L satisfying the stated conditions of the lemma. So, we can consider groups of order rN 
where r divides N, and check if one such group has ordinary irreps with dimensions matching one or more 
of the d\. If we find such an L, and if L has a cyclic subgroup K of order r in its center, then G = L/K has 
a projective irrep of dimension d\ (note that L will have a cyclic subgroup of order r in its center iff this 
irrep includes characters equal to d\, uj r d\, to^dx, • • • , u)* d\). If, after searching groups of order rN for 
all r dividing N, we have not found all the required irrep dimensions, then no group of order N can be used 
for the expansion of Lt, and one should move to groups of larger order. On the other hand, if all required 
irrep dimensions are found, then one must check to make sure that for each central extension L\ (providing 
the projective irrep of dimension d\), the group G = L\/K\ is the same for all A. If it is, one more check 
is necessary, that one can choose projective irreps of the required dimensions that all share the same factor 
system for that particular group G. If so, then we may use the group G to expand hi. 

Alternatively, since the lemma guarantees a single L providing all the irreps that share a single factor 
system for a given G, we can simply search for an L of order rN with irrep dimensions matching the desired 
set. In addition, noting that every normalized factor system may be written as powers of uin (r = N), one 
may choose to restrict the search to groups of order ./V 2 since if the required set of projective irreps sharing 
a single factor system does in fact exist, there must be an extension group of this order providing all these 
irreps. 

The simplest example is for N = 4 when we need only one irrep of dimension d\ = 2. Since every group 
of order \G\ < 6 is Abelian, all such groups have only one-dimensional ordinary irreps. However, setting 
r = 4 so that uj r — i and looking at groups of order rN — 16, we find the group (L) generated by the Pauli 
matrices having a two-dimensional irrep with characters 2, 2i, —2, and — 2i, which tells us that the center 
of this group is isomorphic to the cyclic group, C4, of order 4. In addition, L/C4 — C2 x C2 = G. A valid 
two-dimensional projective irrep for G consists of the identity along with the Pauli matrices, and may be 
readily obtained from the two-dimensional ordinary irrep of L. Note that one can also construct a projective 
representation using r = 2, the central extension then being Q8 or DA (either will work). 
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